KNOTS IN INTERACTION 
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We study the geometry of interacting knotted solitons. The interaction is local and ad- 
vances either as a three-body or as a four-body process, depending on the relative orienta- 
tion and a degeneracy of the solitons involved. The splitting and adjoining is governed by a 
four-point vertex in combination with duality transformations. The total linking number is 
preserved during the interaction. It receives contributions both from the twist and the writhe, 
which are variable. Therefore solitons can twine and coil and links can be formed. 
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Knotlike configurations appear in a variety of physical, chemical and biological sce- 
narios. Examples include Early Universe cosmology, the structure of elementary parti- 
cles, magnetic materials, turbulent fluid dynamics, polymer folding, and DNA replica- 
tion, transcription and recombination. Whenever knots occur, it becomes important to 
understand how they twist, coil, split and adjoin. 

Here we outline a geometrical first principles approach for understanding knots in 
interaction. For definiteness we shall consider Hamiltonian field theories where knots 
appear as solitons i.e. as stable finite energy solutions to the pertinent nonlinear equa- 
tions of motion PJ-[|[]. Even though our conclusions are independent of any detailed 
stucture of the equations of motion, for concreteness we specify to a three component 
vector field n(x) with unit length n • n = 1 ||. We note that n appears as an order 
parameter in a variety of applications. It can for example describe magnetization in a 
ferromagnet, velocity distribution n = v/|v| in fluids, the director field in a nematic 
liquid crystal, and quantum fluctuations of a Higgs field around a spontaneously broken 
vacuum in a grand unified field theory. 

A soliton configuration is localized, and n goes to a constant vector n(x) — * no at 
spatial infinity. Therefore it is a map from the compactified three-space R 3 ~ S 3 to the 
internal unit two-sphere Sg. Such maps fall into nontrivial homotopy classes ^(S ) — Z 
labelled by the Hopf invariant || which can be evaluated by noting that the pre-image 
in R 3 of a point in 5*| is an embedding of a circle S 1 G -R 3 . This embedding is a knot 
in the mathematical sense || , and the Hopf invariant is the linking number of any two 
such knots. We select the asymptotic vector n to point down, along the negative z-axis 
in three-space. At the core (center) of the knotted soliton n then points up, and the 
pre-image of the core is a mathematical knot x(cr) parametrized by a G [0, 27r] so that 
x((t + 2tt) =x(a). 

We consider a short strand of a knotted soliton, a tubular distribution of energy 
density that surrounds the core and does not overlap with any other strand. We inspect 
it with a planar cross sectional disk T> T G -R 3 that cuts the strand at right angle to its 
core at x(er). The center of T> T coincides with the core, and its boundary is a circle 
around the core where n — > no = — e z , see figure 1. Note that topologically each T> T is 
a (Riemann) two-sphere S%, and n has the structure of a map S 1 x S% —>■ 5| that can 
be locally parametrized by 
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Here r is a modular coordinate for the self-linking along the core, it is proportional to 
a. The integer k measures the net number of twists in our strand when we proceed 
along the core. The coordinates tp and 9 are defined at each a, and cover the disk T> T . 
The angle tp g [0, 2ir] is the azimuth around the center of T> T and 6 G [0, it] is the polar 
angle with 9 = at the core and 9 = n at the boundary of T> T . The integer I labels 



the ^(S 2 ) — Z homotopy class of n as a map from the disk T> T ~ S% to the internal 
two-sphere S|. 

Our interpretations become more apparent when we substitute (0) in the integral 
representation of the Hopf invariant Q, 

k ■ I r 
Q H = — - / sin 9d9 A dtp A dr (2) 

Notice that a knotted soliton has either a right-handed or a left-handed orientation 
which is determined by d9 A dip A dr. Furthermore, for any regular Hamiltonian the 
Hopf invariant is conserved even when we have interactions where strands split and 
adjoin. 

We combine n\ and n<i in (|I|) into a two component vector field w = {yi\,n-i). The 
core 9 = is a fixed point of w. Therefore we can define an index Indtp(w) that counts 
how many times w rotates around the core when we go once around it along some curve 
on the disk T> T in the right-handed direction. This index is a topological invariant which 
is independent both of the curve on T> T and of the point a along the core, and from ([!]) 
we find Ind<p(w) = I. This index is also additive, if we evaluate it along a curve that 
surrounds N fixed points of w the result coincides with the sum of the N indices for the 
individual fixed points. 

The third component n 3 in (|l|) is a continuous function on T> T with values between 
[— 1, +1]. It has a maximum at the core 6 = and its minimum occurs on the boundary of 
the disk T> T where 9 — > n. In general the maximum at 6 = is a degenerate critical point, 
and for our knotted solitons we expect that its degeneracy has a definite multiplicity 
given by the absolute value of the index Ind^w), 

deg(n 3 )\ e=0 = |Jnd v (w)| ~ |/| (3) 

In particular, when Ind v (w) = ±1 the core is non-degenerate. In the model studied in 
P]]-[f| we can verify (|3]) explicitly: If p measures distance from the core in local polar 
coordinates ($, p) on T> T1 we find from the equations of motion that the core is truly a 
degenerate critical point of n 3 on D T with /-fold multiplicity, 

nsOM P = l-O(p^) (4) 

This result is related to analyticity near 6 = 0: The disks T> T have a natural complex 
structure sin#<i# A dtp ~ dz A dz suggesting that (|3p should be valid whenever the 
equations of motion are elliptic. Indeed, in our physical applications (^[) makes sense. It 
means that a knotted soliton with Z-fold degeneracy is a bound state of / non-degenerate 
ones. 

Since 9 = is a critical point of n 3 , we have 
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(9 = 0) = (5) 



and since the core is a curve in i? 3 , the 3x3 symmetric matrix 

d 2 n 3 



dxidxj 
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has one vanishing eigenvalue Ai = 0, the corresponding eigenvector is tangent to the 
core. Since 9 = is a maximum, the two additional eigenvalues of (|6]) are non-positive 
A2,3(t) < at the core. If both are non- vanishing the core is a non-degenerate critical 
point of n 3 on T> T . But if the core is a degenerate critical point then either A2 or A3 (or 
both) should vanish. 

We now proceed to apply these considerations to geometrically describe the inter- 
actions between strands of knotted solitons. We first inspect the situation with two 
initial strands that both have a non-degenerate core. We then explain how the results 
generalize when strands have degenerate cores. 

When both initial strands have a non-degenerate core, the corresponding eigenval- 
ues A2,3(t) are non- vanishing. Therefore a nontrivial interaction such as splitting and 
adjoining can not occur unless at least one of the A2,3(t) vanishes. At that point the 
two cores must coincide. 

We first assume that the two initial strands have a parallel relative orientation. From 
(0) we then conclude that when the cores coincide 9 = becomes a doubly degenerate 
critical point of n 3 with Ind v (w) = ±2, depending on whether the two strands are right- 
handed or left-handed. The interaction is pointwise with two outgoing non-degenerate 
strands, and it can be described by a four-point vertex. 

To visualize the vertex we draw plane projections of the strands. For this we trace 
the three gradient vectors dO, dip and dr along the cores. The gradient dr is (co)tangent 
to a core. In a plane projection we describe it by drawing an arrow along the direction 
of the core. For each r the two additional gradients d9 and dip span the cross sectional 
planes T> T . We describe them by drawing two oriented lines in the vicinity of the core. 

We recall that the self-linking number Lk of a knotted soliton can be computed from 
its plane projection using 

Lk = Tw + Wr (7) 

where Tw is the twist and Wr is the writhe in the plane projection. During an interaction 
neither Tw nor Wr is in general separately conserved. Only their sum is since it coincides 
with the Hopf invariant which is conserved. 

To picture the interaction vertex, we first employ continuity to translate all twists 
and writhes in the initial strands away from the vicinity of the interaction point. The 
plane projections of the initial configurations then become well-groomed lines with either 
antiparallel (A) or parallel (P) alignment. For example, in the case of two (right-handed) 
planar strands with antiparallel (A) alignment, we have the four-point interaction vertex 
Va in figure 2. Here both the twist and the writhe are separately conserved. In the case 
of two (right-handed) planar strands with parallel (P) alignment we have the four-point 



interaction vertex Vp in figure 3. Here neither the twist nor the writhe are conserved but 
their sum is, since the interaction preserves the Hopf invariant. Note that a difference in 
the final twist and writhe between the two vertices has important physical consequences, 
it suggests that the strands will twine, coil and link. Indeed, by repeating Vp twice as 
in figure 4, we either continue coiling or form a link, depending on the global geometry. 
This is not equally obvious in the case Va- However, we note that in R 3 a writhe can 
always be continuously deformed into a twist. Furthermore, two non-coplanar strands 
can exhibit an antiparallel alignment in one plane projection, but parallel alignment 
in another plane projection. Consequently we expect that the two vertices Vp and Va 
should be related. For this we first vertically flip the plane projection of one of the 
two initial strands in figure 3. This yields two strands with antiparallel alignment in 
the plane projection. We then implement an interaction described by the vertex Va in 
figure 2, followed by another vertical flip. The result coincides with that obtained by 
implementing the vertex Vp alone. This means that the two vertices Va and Vp are dual 
to each other, with vertical flip defining the duality transformation. 

Depending on the dynamical details of the Hamiltonian, besides splitting and ad- 
joining the two initial strands may also combine into a single strand but with a doubly 
degenerate core. This can be described using the topologically invariant index Ind 9 {w). 
As an example we consider two initial right-handed strands both with Ind 9 (w) = +1 
(figure 5). Since the index is additive, we conclude that for a curve which encircles both 
cores we have Ind v (w) = +2. This is also the index of the final, doubly degenerate 
strand. The interaction is a process where two homotopically nontrivial maps S* —*■ S 1 
each with a 7r 1 (S' 1 ) winding number +1 combine into a single map with 7r 1 (S' 1 ) winding 
number +2. Dynamically, this leads to the formation of a doubly degenerate knotted 
soliton, a bound state, in a manner which is consistent with the conservation of the Hopf 
invariant. We visualize the interaction as a three-body process in figure 5. 

When the two initial strands have an opposite orientation, splitting and adjoining can 
not occur since two strands with different orientation should not connect. Instead the 
two strands can annihilate each other. As in figure 5, this can also be described using 
the topologically invariant index Ind 9 {w). For example, if one of the initial strands 
corresponds to Ind v (w) = +1 and the other strand to Ind^w) — — 1, since the index 
is additive we conclude that for a curve that encircles both cores we have Ind^iw) = 0. 
As in figure 5, the interaction is a process where two homotopically nontrivial maps 
S* — ► S* 1 combine. But now one of the initial maps has 7ii(S r ) winding number +1 
and the other — 1, and the combined map is homotopically trivial. This corresponds 
to an annihilation between the two strands, that proceeds dynamically in a manner 
which is consistent with the conservation of the Hopf invariant. The annihilation can be 
visualized as a three-body process, much like in figure 5. 

In the general case where the initial strands are degenerate the interaction can pro- 
ceed in a multitude of different fashions, details depending on the degrees of degeneracy 



and the relative orientations of the interacting strands. But according to (|3]) degenerate 
strands can always be viewed as bound states of non-degenerate strands. Thus there is 
no need to discuss their interactions separately. We can always describe their interac- 
tions using various combinations of the interactions that occur between non-degenerate 
strands, in a rather obvious manner. 

Finally, the interaction processes that we have outlined here all conserve the Hopf 
invariant. This follows from the condition that the order parameter n is a unit vector. 
But in a number of applications one may wish to relax the constraint n • n = 1 so that 
|n| can vanish at some points in R 3 . One way to achieve this is to replace the constraint 
by a term A(n- n — l) 2 in the Hamiltonian, with A a coupling constant. When A —>■ oo we 
recover n-n = 1. But for finite (large) coupling we expect |n| to fluctuate in the vicinity 
of |n| « 1 even though we can also have |n(x*)| — >• at some points. This generalizes the 
models discussed here in a manner that should be of interest in a variety of applications. 
For example in quantum field theory it allows us to recover renormalizability, while in 
molecular biology it enables us to account for effects of entzymes in DNA. Obviously 
study of interaction processes with a finite A should be of wide interest. 

We have found that depending on the relative orientation and degeneracy of the 
solitons involved, an actual interaction can involve various combinations of splitting and 
adjoining, twisting, coiling and linking. It can lead to the formation of bound states and 
annihilation between strands. Our description of the interaction processes is quite model 
independent and based entirely on elementary conceps of continuity and differentiability. 
Consequently we hope that our results can form a basis for detailed investigations of 
interacting knots. Indeed, we expect that numerical studies in the model discussed in 
M-[§j (and its generalization to the case where |n| can vanish) provide prolific tests for 
various aspects of knot interaction. 
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Figure Caption 



Figure 1: A cross-sectional disk T> T cuts a strand at a right angle to its core. The 
center of the disk coincides with the core where n points up, and its boundary is a circle 
where n points down. 

Figure 2: The plane projected interaction vertex Va between two right-handed, 
mutually antiparallel strands. 

Figure 3: The plane projected interaction vertex Vp between two right-handed, 
mutually parallel strands. 

Figure 4: If repeated twice, the vertex Vp can lead to further coiling or to the 
formation of a link. 

Figure 5: Two right-handed non-degenerate strands can combine into a single 
doubly-degenerate strand. The process leads to a doubly degenerate knotted soliton. 
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